Freezing transition of the vortex liquid in anisotropic superconductors 
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We study the solid-liquid transition of a model of pancake vortices in laminar superconductors using 
a density functional theory of freezing. The physical properties of the system along the melting line 
are discussed in detail. We show that there is a very good agreement with experimental data in 
the shape and position of the first order transition in the phase diagram and in the magnitude 
and temperature dependence of the magnetic induction jump at the transition. We analyze the 
validity of the Lindemann melting criterion and the Hansen- Ver let freezing criterion. Both criteria 
are shown to be good to predict the phase diagram in the region where a first order phase transition 
is experimentally observed. 
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I. INTRODUCTION 

The thermodynamics and dynamical properties of vor- 
tex matter in high temperature superconductors have re- 
ceived the ijiterest of a large experimental and theoretical 
communityld. One. of the most salient findings, already 
predicted in 1985a, is the first order melting transition 
of the vortex lattice. Large thermal fluctuation in a high 
Tc material induce the melting transition first observed 
in transport measurementsQ and confirmed by a number 
of equilibrium propssties like direct nieasurements of the 
magnetization jumpQ and latent heatS. 

It is now well accepted that in clean Bi2Sr2CaCu208 
(BSCCO) and untwined YBaaCuaOr (YBCO) samples 
the melting of the vortex lattice for an external magnetic 
field parallel to the c-axis of the samples is first order 
and that the superconducting coherence is lost simulta- 
neously in all directions at the melting temperature. De- 
spite of the great amount of work devoted to the subject, 
the nature of the melting transition of the vortex lattice 
is still controversial. 

The first theories for melting of the vortex lattice were 
based on the Lindemann criterionfl and describe a situ- 
ation in which the vortex lattice melts into a liquid of 
entangled vortex lines. However the behavior of the en- 
tropy change at the melting transition suggest that at the 
first order transition a simultaneous loss of the triangu- 
lar crystal structure and a decoupling of the planes take 
place. In this last scenario, the liquid would be a liquid of 
pancake vortices. Many authors have analyzed the exper- 
imental results obtaining good fittings with this picture. 
The precise nature of the transition and the structure 
of the liquid phase at the melting line is still an open 
problem. 

More recently, some striking regularities -weak 
universalities- observed in the transport properties at 
the melting tegwerature were interpreted in terms of the 
Hansen- Verlell3,B freezing criterion. This criterion, that 
has been proved to work very well in the case of classi- 
cal liquids, states that the liquid freezes when the first 



peak of the structure factor S{k) reaches a critical value. 
However, for the case of vortex matter, the validity of 
this criterion has not been analyzed. 

In this context it is interesting to revisit the density 
functional theory for vortices in the highly anisotropic 
high Tc superconductors. This theory describes inter- 
acting pancake vortices in layered superconductors. The 
interactions are of electromagnetic nature, Josephson in- 
teractions involving charge transfer are neglected, and 
consequently it is a good starting point to describe 
highly anisotropic systems. However, as we discuss be- 
low, it may be also appropriate to describe even the less 
anisotropic high Tc cuprates. The theory allows to calcu- 
late the melting temperature, thermodynamic properties 
at the melting line as well as the validity of the phe- 
nomenological criteria for melting and freezing. 

In what follows we present the model for the inter- 
vortex interaction and the method for the calculus of 
the melting line. We present the results for the phase 
diagram and compare them with experimental data. We 
analyze the validity of the Lindemann and Hansen and 
Verlet criteria for melting and freezing. 



II. DESCRIPTION OF THE MODEL 

Our starting point is the functional density theory for 
pancake vortices in a layered superconductor in a mag- 
netic field perpendicular to the layers. In a layered sys- 
tem in the limit of infinite effective mass perpendicular 
to the layers (zero Josephson coupling) the pancake vor- 
tices can be treated as point like classical particles re- 
stricted to move in the planes. The vortex-vortex inter- 
action is given by a three-dimensional (3D) anisotropic 
pair potentiafl which in Fourier space is given by: 
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here, k = (kj^, kz) is the wave vector, A is the planar Lon- 
don penetration depth and T = /3d$Q/47rA is a dimen- 
sionless strength parameter where $o is the flux quan- 
tum, d is the distance between planes and /3 = (ksT)"^ 
is the inverse temperature. 

Following the, pioneering work of Ramakrishnan and 
YussoufF (RY)E3 |aH|i|the extensions made to describe 
the vortex mattcrtil,ll3, the functional describing the free 
energy difference between the solid and liquid phases is 
given by: 
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here, r = {r^^nd) where rj_ is the in-plane coordi- 
nate and n is an integer (the plane index). The quan- 
tity Pn(i"_L) =< ^i5{v±_ — r"-) > is the vortex den- 
sity with r" j the 2D coordinate of the i-th. particle at 
plane n, pi is the mean areal density of the liquid and 
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I) is the direct correlation function in 



the liquid phase. The first integral in the right hand 
side of Eq. ^ describes the non-interacting contribution 
to the free energy, the second integral incorporates the 
effect of the interactions up to second order in the differ- 
ence (p„(r^) — p<').To evaluate this energy difference it is 
necessary to calculate the direct correlation function of 
the liquid defined below. 

The pair distribution function is defined as: 
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where: 



here pn,ra'(r±,r^) is the probability density of finding 
two particles at r —{r±,nd) and r'^{r'j_, n'd) respec- 
tively. Due to the symmetry of the liquid g(r,r') = 

9\n-n'\i\r±-T'j_\). 

The direct correlation function (pr/rj.) is defined by 
means of Ornstein-Zernike equationE3: 

hn{rj_) = c„(r_L) + Pi"^ j d^r'^c\„-n'\{\ri. - ?'j_|)^n'(f'j_)- 

n' 

(5) 

where hn{r±^) = 9n{'>'\_) — 1 is the pair correlation 
function. 

This integral equation for the direct correlation func- 
tion c„(r_L) can be solved in the hipernetted chain ap- 
proximation, and inserting the solution in Eq.(|^), the 
free energy difference between the solid and liquid phase 



is given only as a functional of the density Pn(r_L ). At 
this point, we can expand the density in a set of trial 
reciprocal-lattice vectors: 
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The quantities rj and , may be viewed as dimension- 
less order parameters, rj is the fractional density differ- 
ence between the liquid and the solid phases and the pc 
are the Fourier components of the density in the solid 
phase. Minimizing the free energy with respect to 
these parameters, we obtain the following equation: 

1 + + ^ pcje''^'' = cxp(p£Co77 + pi^ pGCGe''^'"), 

G^O G^O 

(7) 

where cg is the Fourier transform of the direct correlation 
function and cq = cg=o- 

Due to the small compressibility of the system the frac- 
tional density change is small (jj <ti 1), and as a first 
approximation can be neglected in the first term of (|^). 
Note that the product cqT] is not necessarily small and 
the first term in the exponential of the right hand side of 
(0) has to be retained. Transforming Fourier Eq. (||) we 
obtain after some algebra a set of coupled equations for 
the pg- 
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here A^, is the area of the in-plane unit cell. Given the 
Fourier components cg of the direct correlation function, 
these equations can be solved for {pg}- Before present- 
ing the numerical results, the following points deserve a 
comment: 

i) A uniform liquid is always a solution of (^) {pG = 0}. 
As the Cg increase, non-uniform solutions may suddenly 
appear (with pG ''^ 0.5) and correspond to the solid phase 
with a spatially periodic density. 

ii) At the transition the free energy difference Afi de- 
fined in ^ is equal to zero. To calculate the transition 
line in the B — T phase diagram, we have to calculate the 
direct correlation function as a function of the magnetic 
induction and temperature and use it as an input in the 
RY density functional theory. The transition is obtained 
when a set of order parameters {pG 7^ 0} given by (||) 
satisfy the condition Ail — 0. 

iii) As discussed in previous works, the numerical cal- 
culation can be simplified by noting that in the liq- 
uid phase, the Fourier transform cg of the direct cor- 
relation function rapidly decay with increasing |G|, it 
is then a good approximation to truncate the series 
Eg^-^o PGCce**^ '' in (^) and take only a few terms with 
Cg 0. The number of terms retained defines the num- 
ber of self consistent order parameters {pg} included in 
the theory. In what follows we consider only two order 
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parameters corresponding to the first two peaks of the 
in-plane direct-eorrelation function as first suggested by 
RamakrishnantJ for a 2D system. 

iv) With the resulting set of order parameters that 
characterize the solid phase at the melting point in the 
limit T] ^ I, the fractional density change can be esti- 
mated. However in the RY approximation used here the 
results for rj are not accurate, in particular it is ease to 
show that 7/ > in contradiction to what is expected for 
the case of vortices. Rather than improving the estimate 
of 77 by including higher order terms in the functional Af2 
of Eq. ^ -which are given by not available three particle 
correlation functions- we estimate the entropy change at 
the transition and calculate t] by means of the Clausius- 
Clapeyron relation. The higher order corrections to Afl 
are not strongly coupled to the order parameters {pc} 
and the present approximation is known to give accurate 
results when the crystal structure of the solid phase is 
given. 



III. CORRELATION FUNCTIONS 

In this section we present results for the pair correla- 
tions in the liquid phase, we use the hipernetted chain 
approximation (HNC) which is well described in the lit- 
erature. For the sake of completeness we only give here 
the basic formulas of classical liquids theory relevant for 
our case. 

The pair correlation function can be written in the 
form: 



^ g[9«('-i)-l-/5V„(r_L)-c„(r_L)+S„(r_L)] 



(9) 



where Bn{r±) is known as the bridge function and its 
closed expression is not known for an arbitrary poten- 
tial. In the HNC approximation we have Bn{r±) = 
which is an excellent approximation for ktag ranged soft 
potentials as the one we are dealing withEJ. In the HNC 
approximation Eq. can be rewritten in the following 
form: 



'PV„ir^)+Y„{r^)-l-Y„{r^)] 
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where y„(rj^) = hn{r±) — Cnir±) . This equation together 
with the Ornstein Zernike relation (Eq.(|^)) are the ba- 
sic equations that allows to evaluate the pair correlation 
functions. The asymptotic behavior of the direct corre- 
lation function as rj^ — > 00 and n — > 00 is known to be 
Cn{r±) = —Vn{r±_)/kBT. For the numerical calculation 
it is then convenient to divide the correlation functions 
in short and long ranged contributions according to: 



with 
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and 



with 
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Due to the highly anisotropic nature of the pair po- 
tential, as a first approximation we take the short ranged 
contribution different from zero only for the in-plane cor- 
relations: 



(15) 



and the problem reduces to calculate a single function 
(cos(''-l) ) given by the following set of integral equations: 

cosir±) = expiYosir^)) - 1 - Yosir^), (16) 
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cos(fc±) - PVjk) 
1 - Pi[cosik±) - PV{k)] 



(17) 



the last equation is the Fourier transformed Ornztein- 
Zernike equation. This expression can be further simpli- 
fied by an analytic integration due to the particular form 
of the pair potential. 

Finally the out of plane short ranged correlation func- 
tions in a second order approximation are given by the 
following equation: 
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In Fig. we present results for the pair distribution 
function where the lengths are in units of the characteris- 
tic distance a given by na^pi = 1. The first peak of the in 
plane pair distribution function occurs at a finite value of 
the which is given by the mean inter-particle distance 
(^ 2a). For the out of plane distribution function the 
maximum occurs at = which shows the tendency of 
pancake vortices to form stacks. In Fig. || the behavior 
of h(r± = 0, n ^ 0) is presented and shows an exponen- 
tial decay as a function of n. This behavior allows as to 
define a correlation length A through the following rela- 
tion: h{r± — 0,n) (X exp(— dn/A). In Fig. ^ we present 
the behavior of A as a function of the magnetic inductioa. 
B, that can be fitted with A (x 1/Vb. Other authorsllll 
had an excellent fit of their data in a smaller range of 
fields with a function of the form A = Ai + A2/B. 

It is important to note that the characteristic length 
A describes the long distance behavior of the correla- 
tion functions. The short distance behavior is given by 
h{rj_ = 0, n) with small n which in general decays very 
fast. As can be seen in Fig.pl h{r± = 0,n = 1) is very 
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small (~ 10~^), and although its value depends on field 
and temperature, in the liquid phase, it always remains 
much smaller than 1. 

As the temperature decreases the system develops 
stronger correlations. This is clearly observed in the 
structure factor, defined by: 



,2 „2 



_ -\-kz nd) 



(19) 

In Fig. ^ the structure factor is shown for two different 
temperatures, the lower the temperature, the higher the 
first peak in 5'(k;). 

In the next section we use these results to calculate the 
melting line. 
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If the characteristic length / is taken as the distance be- 
tween planes which should be a good approximation at 
least in the high field regime we obtain: 



ksTm 
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(22) 



If the inter-plane correlations were relevant at the melt- 
ing point, the characteristic length / should be obtained 
minimizing the free energy (21). In this case the melting 
line is given by: 



(23) 



IV. PHASE DIAGRAM 

We used the correlation functions calculated above as 
an input to the RY density-functional theory of freezing. 
In Fig. [s] we show the results obtained for the phase 
diagram with parameters corresponding to BSCCO but 
without temperature dependence: A(T) — A(0) = ISOOA, 
d = IsA and Tc — 90K. To understand these results we 
resort to a simple model obtained by mapping the vortex 
system onto a 2-D boson system. In this .iBodel, the free 
energy of a system of vortices is given bytS: 



dz 



(20) 



where rj(z) determines the position and shape of the i-th. 
vortex, g is the elastic energy and the last term describes 
the interaction between vortices, which is taken to be 
logarithmic, rij(z) is the distance between vortices. In 
a solid phase, of lattice constant ao, we can use an har- 
monic approximation for the pair interaction potential. 
Assuming a characteristic length I which defines the cor- 
relation along the field direction, the free energy after a 
thermodynamic average can be approximated by: 



F' = F/N 



< > 



(21) 



where F' is the mean free energy per characteristic vor- 
tex segment of length I, N is the number of independent 
segments in the sample, < > is the mean square dis- 
placement of the vortices and k cx eo/oo curva- 
ture of the logarithmic pair potential at the mean inter- 
particle distance. Using the Lindemann melting criterion 
< >= u|aQ, with ul ~ 0.1 — 0.3 along the transi- 
tion line and the equipartition theorem that states that 
F' ~ fc^T, the expression above reduces to: 



The density functional data for the melting line were 
fitted with the following expression Tm = + Ai/Bm 
according to ([2^). This type of behavior is expected for 
a highly anisotropic system at high fields, when the liq- 
uid state at the melting temperature presents very short 
correlations along the c axis. Although the fit is less ac- 
curate for small fields we never obtain a behavior of the 
type given by (|23|). This suggests that in our field and 
temperature range, the system simultaneously melts and 
the planes decouple to form a liquid of pancake vortices. 
The less accurate fit for the low field data could indicate 
that we are entering a cross-over to a line melting regime. 
In the limit of high fields the melting temperature goes 
asymptotically to the 2D melting temperature T^. This 
is due to the fact that at high fields the inter-plane cor- 
relations become irrelevant and the system behaves as a 
collection of independent planes. The 2D melting tem- 
perature, obtained with a functional density theory in 
the RY approximation, is lower than expected from a 
Kosterlitz-Thouless dislocation unbinding theory. This 
is presumably due to the fact that the HNC approxima- 
tion underestimates the correlalions in the 2D case as 
was discussed by other authoraHl. 

In order to compare our results with the experimen- 
tal data, we includpd the temperature dependence of the 
penetration lengthM: A(r) = A(0)/[1 - (T/Tc)'^]^^^ . In 
Fig. H we present our results with A(0) — 1500A to- 
gether with the experimental results of Zeldov et al. In 
the experimental data the first order transition ends at a 
critical point presumably due to point-like pinning cen- 
ters, an effect not included in this functional density the- 
ory. In order to make this comparison, A(0) was used as 
a fitting parameter. A more realistic value of A(0) for 
BSCCO is - 2000A. If the A(0) is increased, the melt- 
ing line moves to lower fields due to the decrease of the 
intensity of the pancake-pancake interactions, this could 
be compensated by the inclusion of a small Josephson 
coupling between planes in a more realistic treatment. 
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V. ENTROPY AND MAGNETIC INDUCTION 
JUMP AT Tm 

The entropy jump can be calculated as the tempera- 
ture derivative of the free energy difference: 



A5 



(24) 



Using Ecs. H and ^ the free energy difference of Ec| 
can be written as: 



In 
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(25) 
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To evaluate the entropy jump we must calculate the 
temperature derivatives of the Fourier transforms of the 
liquid direct correlation functions cg at constant H and 
V , that are not known. Following referenced, this tem- 
perature derivatives can be estimated on the solid phase 
using the fluctuation-dissipation theorem in the classical 
limit. This procedure gives for the structure factor S{G) 
a linear dependance on temperature. Using this result 
we obtain the following expression for the entropy jump: 



As/ZcB = 3[p?(l-ci)+p2(i_c2)]. 



(26) 



in our two order parameter calculation. 

In Fig. ^ we present the density functional results for 
the entropy change. Its average value ~ Q.iks is con- 
sistent with experimental datau. With the entropy jump 
result and the Clausius-Clapeyron relation we can calcu- 
late the magnetic induction jump at the transition line: 
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/ d<i>o dB„, 
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where we have made the following approximation: 



dT 



dB„i 
~dT 



(27) 



(28) 



The results for the relative magnetic induction jump are 
plotted in Fig. |^ with the Zeldov et al. experimental data 
digitalized. The parameters used to obtain this data were 
the same as those used to calculate the phase diagram of 
Fig. p. There is a very good agreement in the magnitude 
and temperature dependance with this direct experimen- 
tal measurement. 



VI. MELTING AND FREEZING CRITERIA 



melts when the square root of the mean square displace- 
ment of the particles is an order of magnitude smaller 
than the lattice parameter i.e. < > /ag = with 
Ul 0.1. This criteria has been widely used to estimate 
the phase diagram and different properties of the system 
at the melting line. In general it is accepted that for the 
melting of vortex mater the Lindemann parameter is 
larger than for other classical 3D solids and its value is 
usually taken in the range 0.2 — 0.3. 

We estimated the mean square displacement as: 



1 



dWrl^^. 
Pi 



(29) 



the density at the melting line is obtained using (^ and 
the results for the Lindemann parameter are shown in 
Fig. ^. For low fields, numerical results show that ul 
~ 0.2 in good agreement with previous estimates, as the 
field is increased (the melting temperature decreases) the 
Lindemann parameter increases. This is due to the fact 
that, as discussed above, in the present model as the 
field increases the interplane coupling decreases and the 
transition becomes essentially 2D. For a 2D transition, 
the Lindemann parameter is larger than for a 3D one. 

The Hansen- Verlet freezing criteria is the counterpart 
of the Lindemann criteria and states that a liquid freezes 
when the first peak of the structure factor reaches the 
value 2.85 in a 3D classical system and of the order of 
5.5 for a 2D plasmallJ. Having evaluated the structure 
factor at the melting transition within the functional den- 
sity theory, we can show the validity of this criteria for 
the case of freezing of the vortex liquid. The results for 
'S'(fcmax) are shown in Fig. ^ and as in the case of the Lin- 
demann criteria the Hansen- Verlet criteria is valid for low 
fields with a value 5'(fcmax) ^ 6.5, much larger than that 
obtained for isotropic 3D system. As the melting tem- 
perature decreases >S'(fcinax) towards its two dimensional 
value. In our case the value of 5'(fcinax) at small fields 
is larger than the 2D value due to the following reason: 
the idea behind the Hansen- Verlet criteria is that when 
the liquid develops correlations, that overcame some crit- 
ical value, it freezes. At high temperatures the first peak 
in S{k) dominates and is a good measure of the corre- 
lations, as the temperature decreases, the liquid devel- 
ops stronger correlations for large G and both the first 
and second peaks in S{k) became relevant at the freezing 
point. Consequently, it is not necessary for 5(fcniax) to 
become too large to indicate that the liquid has developed 
critical correlations. This can be viewed as a relative in- 
crease of the order parameter p2 with respect to piat the 
transition as the melting temperature decreases. The re- 
sults indicate that for low fields the ratio pi/p2 is larger 
than 4, while at high fields is of the order of 2. 



In this section we analyze the Lindemann melting crite- 
rion and the Hansen- Verlet freezing criterion for the vor- 
tex system. The Lindemann criteria states that the solid 
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VII. SUMMARY AND DISCUSSION 

In this work we study the melting transition for the 
vortex system in layered superconductors using the RY 
functional density theory. Using reasonable set of param- 
eters for BSCCO we obtained qualitative and quantita- 
tive agreement with the available experimental data. 

The functional density theory allows an estimation 
of the entropy jump at the transition and using the 
Clausius-Clapeyron relation we estimated the jump in 
magnetic induction. The same parameters used to fit 
the phase diagram, give results for the magnetic induc- 
tion jump that are in quantitative agreement with the 
experimental results. Finally we analyzed the validity of 
the melting and freezing criteria and showed that in the 
region of interest -the low field region where the first or- 
der transitions occurs- the Lindemann and the Hansen- 
Verlet criteria are valid within a error of 10% which is 
typical for this type of criteria. 

The model predicts the melting of the vortex lattice 
into a liquid of almost decoupled pancake vortices even 
in the low field region. In this sense we may interpret 
our results as a suaultaneous decoupling and melting of 
the vortex latticeEj rather than a melting into a liquid 
of vortex lines. This picture may be appropriate for a 
highly anisotropic system like. BSCCO. However, in a 
recent paper Kitazawa et al.E^ successfully used a scal- 
ing suitable for simultaneous decoupling and melting not 
only for BSCCO but also for less anisotropic systems like 
YBCO. This scaling predicts that the melting line is given 
by: 



B„ 



. J- in. 



(30) 



where the prefactor depends on the anisotropy and inter- 
plane separation of the system. Both the experimental 
results of Zeldov et al. for BSCCO and our numerical re- 
sults for low fields are in good agreement with this expres- 
sion as shown in Fig. ^ In the figure we present results 
corresponding to a zero temperature penetration depth 
A(0) = 1500A which gives very good agreement with the 
results of Zeldov et al. for the low field region where 
the first order transition is experimentally observed. As 
discussed in section IV, a more realistic value for A(0) is 
2000 A, which gives a melting line laying below the ex- 
perimental points. However, in any real system there is 
some Josephson coupling between the planes that should 
compensate the lowering of the electromagnetic coupling 
due to the increase of A(0). In this sense, we take A(0) 
as an effective parameter that measures the strength of 
the coupling between planes. A change in A(0) does not 
change the temperature dependance of B„i for low fields, 
but only the coefficient C in Eq.(|30|). 

If experimental results confirm the fact that even for 
less anisotropic system, the first order melting transition 
correspond to a simultaneous decoupling and melting of 



the vortex lattice, the present theory could be used to de- 
scribe these systems provided that A(0) or the interplane 
distance are chosen as effective parameters. In particular 
the validity of the melting and freezing criteria obtained 
for the highly anisotropic laminar BSCCO could be ex- 
tended to other less anisotropic systems like YBCO. 



^ G. Blatter et al, Phys. Rev. B, 54, 72 (1996); R. Cubitt et 

al, Nature (London)., 365, 407 (1993). 
^ E. Brezin, D.R. Nelson, A. Thiaville, Phys. Rev. B , 31, 

7124 (1985). 

^ H. Safar et al, Phys. Rev. Lett. 69, 3370 (1992). 

■'E. Zeldov, D. Majer, M. Konczykowsky, V. G. Geshken- 

bein, V.M. Vinokur, H.Shtrikman, Nature, 375, 373 

(1995). 

^ A. Schilling et al. Nature 382, 791 (1996). 

A. Lindemann, Z. Phys. 11, 609 (1910). 

J. P. Hansen, L. Verlet, Phys. Rev. 184, 151 (1969). 
^ A.G. Rojo, A.K. Sood, C.A. Balseiro, Solid State Comm. 

,107, 401 (1998). 
^ Feigel'man, et al, Physica C 167, 177 (1990). 
^° T.V. Ramakrishnan, M.Yussouff, Phys. Rev B 19, 2775 

(1979). 

" G. I. Menon et al, Phys. Rev. B 54, 16192 (1996). 

S. Sengupta, C. Dasgupta, H.R. Krishnamurthy, G.L 

Menon, Phys. Rev. Lett. 67, 3444 (1991). 
^■^ J. P. Hansen, I.R. McDonald, Theory of Simple Liquids, 

Academic Press, London (1986). 
" T. V. Ramakrishnan, Phys. Rev. Lett. ,48, 541 (1982). 

G.W. Crabtree D.R. Nelson, Physics Today, 50, 38 (1997). 

S. L. Lee et al, Phys. Rev. Lett. 75, 922 (1995). 

P. Vashishta, R.K. Kalia, Proceedings of the Ninth Midwest 

Solid State Theory Symposium, U.S.A. (1982). 
1* L.I. Glazman and A.E. Koshelev, Phys. Rev. B 43, 2835 

(1991); L.L. Daemen et al, Phys. Rev. Lett. 70, 1167. 

(1993). 

T. Sasagawa. K.Kishio, Y. Shimoyama, K. Kitazawa, Phys. 
Rev. Lett. 80, 4297 (1998). 

Work partially supported by grants CONICET 
4946/96 and ANPCYT PICT 0251/98. One of us 
(C.A.B.) is partially supported by CONICET. 



6 




B (Tesla) 



FIG. 1. HNC results for the pair distribution function 

g{r±,n = 0,1,2,...) of the vortex liquid for B = 3T and 
T = 17K. For n = 0, the first maxima as a function of r± is 
given by the mean particle spacing and for n is at r± = 0. 



FIG. 3. The c-axis correlation length as a function of the 
magnetic induction B with T = 31K. The slope of the linear 
fit in the inset is 1/2. 



0.015 




5 10 15 20 25 30 

n 

FIG. 2. The pair correlation h{r±,n) at rx = as a func- 
tion of n for B = 0.065T and T = 31K. The values can be 
fitted with an exponential decay, as shown in the inset with 
a log- normal scale. 
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FIG. 4. Structure factor for kz — 0. Note the increase in 
the structure as the temperature is lowered. The amplitude 
of the first peak is the quantity that we considered for the 
Hansen- Verlet criterion. 
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FIG. 5. Liquid-Solid transition line in the B-T phase dia- 
gram for the vortex system in BSCCO. The solid line is fit 
with B = Ai/{T - T^), where Ai and are variational 
parameters. 



FIG. 7. Entropy jump at the first order transition calcu- 
lated with the RY functional density of freezing. 




FIG. 6. Phase diagram of BSCCO. The open dots are the 
RY functional density results and the solid dots are direct 
meassurements by Zeldov et al. In the experimental data 
there is a critical point where the first orther transition ceases 
to exist, which is prossumably caused by the presence point 
pinning centers in the sample. 




FIG. 8. Relative magnetic induction jump at the first or- 
der transition. The open dots wore calculated using the Clau- 
sius-Clapeyron relation and the RY functional density theory 
of freezing. It is in very good agreement with the experimen- 
tal results of Zeldov et al. (solid dots) in the magnitude and 
in the temperature dependence. 
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FIG. 9. Lindemann melting parameter and Hansen- Verlet 
freezing parameter as a function of Tm for BSCCO calculated 
in the RY approximation. Both parameters are approximatly 
constant in the region where the first order phase transition 
it is experimentally observed. 




FIG. 10. Bm asa. function of Tc/T - 1. A test for subli- 
mation. Calculations with different A(0) or d give liquid-solid 
transition lines displaced paralel to each other giving the same 
effect of a different anysotropy. 
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